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^L(' Abstract 

«]-> : 

(~| , Direct detection of dark matter (DM) requires an interaction of dark matter particles with 

nucleons. The same interaction can lead to dark matter pair production at a hadron collider, and 

*vj ! with the addition of initial state radiation this may lead to mono-jet signals. Mono-jet searches at 

^ ' the Tevatron can thus place limits on DM direct detection rates. We study these bounds both in the 

C^^ \ case where there is a contact interaction between DM and the standard model and where there is a 

0^ ■ mediator kinematically accessible at the Tevatron. We find that in many cases the Tevatron provides 

the current best limit, particularly for light dark matter, below ~ 5 GeV, a and for spin dependent 
interactions. Non-standard dark matter candidates are also constrained. The introduction of a light 

|0 \ mediator significantly weakens the collider bound. A direct detection discovery that is in apparent 

conflict with mono-jet limits will thus point to a new light state coupling the standard model to 
the dark sector. Mono-jet searches with more luminosity and including the spectrum shape in the 
analysis can improve the constraints on DM-nucleon scattering cross section. 
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1 Introduction 

From astronomical and cosmological observations it is now clear that ^ 25% of the matter-energy 
content of the universe if made up by dark matter (DM). Although DM has so far only been observed 
through its gravitational interactions the quest for a more direct observation of DM is taking place 
simultaneously on many fronts. Indirect searches look for signals of standard model (SM) particle 
production from DM annihilations in our galaxy, direct searches look for interactions of DM with SM 
particles in underground detectors and colliders attempt to produce the DM and measure it. We will 
concentrate here on direct detection and collider searches. 

If dark matter is to be observed in direct detection searches it must couple to quarks or gluons ^ . 
The same couplings lead to direct DM production at hadronic colliders such as the Tevatron, and 
we wish to investigate the connection between the two types of search. We will do so in a model 
independent fashion [1] ; we will assume that the DM is fermionic and that there is some massive state 
whose exchange couples DM to quarks. The mediator may be a SM gauge boson, the Higgs or a new 
particle (if the new particle is very heavy we can describe its effects with an effective contact operator). 
Although the processes that give direct detection and those that give DM production occur through 
s- and t-channel exchange of the same mediator, the regimes probed in the two types of experiment 
are very different. The momentum exchange during a DM-nucleus recoil is ~ 100 MeV whereas at the 
Tevatron the typical momentum exchange is 10 — 100 GeV. This leads to two interesting regimes to 
consider when comparing bounds from the two types of experiments: heavy mediators M > 100 GeV 
and light mediators M < 100 GeV. 

The momentum exchange at direct detection experiments is sufficiently low that for all but the 
lightest mediators below 0(100 MeV), which we do not consider here, the mediator can effectively be 
integrated out and the scattering rate in both regimes scales as, 

CTDD~ffx5g^> (1) 

where, for simplicity, we have ignored form factors and possible momentum and velocity dependence 
in the cross section. Here, g^ and Qq are couplings of the mediator to DM and quarks, n is the reduced 
mass of the DM-nucleon system. 

In contrast the two regimes behave very differently at colliders. Concentrating on direct production 
of a pair of DM particles and an initial state emission of a jet, we estimate the mono-jet + fJx partonic 



'^DAMA is an exception as, unlike other experiments, it does not distinguish between nuclear and electron recoils of 
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production cross section in the two cases to be 
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where as is the QCD couphng and pT is the transverse momentum of the jet, which is typicahy 
~ 0(100) GeV at the Tevatron^. Thus, for the heavy mediator case the (partonic) production cross 
section at the Tevatron, where px ~ 100 GeV, is 0(1000) times larger than the direct detection cross 
section for /i ~ 1 GeV when the DM is heavier than the nucleon mass. The CDF mono-jet search [2] 
analysed r^ 1 fb^ and saw no significant discrepancy from the SM, thus limiting the DM + mono-jet 
production cross section to be smaller than ~ 500 fb. Due to the factor of 1000 mentioned above, 
this will translate to bounds in the neighborhood of 0.5 fb in direct detection experiments, the exact 
bound at direct detection experiments will depend upon the details of the parton density functions 
relating the partonic cross section of (2) to the actual CDF mono-jet bound. 

This is to be compared with direct detection current searches. Null results from experiments such 
as CDMS [3], XEN0N[4, 5] and others, place strong constraints on the cross section of DM to recoil 
from a nucleus, a < 10^^ — 10^^ fb for a 10-100 GeV WIMP scattering elastically through a spin 
independent (SI) interaction. Thus, for this situation it seems that direct detection has greater reach. 
However, due to the threshold to detect a DM recoil in these experiments there is a DM mass below 
which these experiments are no longer sensitive, typically this lower bound is m^ ~ 5 — 10 GeV, there 
is no such threshold in collider searches. 

Furthermore, the DAMA collaboration [6] have observed a signal consistent with DM scattering 

from Nal which is inconsistent with bounds on a standard WIMP from CDMS and other experiments. 

This has motivated the introduction of non-standard DM scenarios that can make these seemingly 

discrepant results consistent. The cross sections necessary to explain DAMA are considerably larger 

than 10~^fb and may allow these scenarios to be probed directly at the Tevatron, due to the increase 

in cross section described above. Another possibility that has been motivated both by DAMA and 

the recent CoGeNT [7] excess is that dark matter is light, below about 10 GeV, and is thus transfers 

small momenta to nuclei giving a signal near threshold. The Tevatron will place a strong bound for 

dark matter particles below 5 GeV. Finally, spin-dependent (SD) WIMP-nucleus scatterings are not 

coherent and therefore are not enhanced by an A^ factor. Typical bounds on a SD WIMP-proton 

scatter from direct detection are ~ 1 fb , and will be severely impacted by the mono-jet bounds 

presented here. 

^Note that (2) is only qualitative in nature. The limits are correct for mediator masses well above and below the pr 
of the jet. 



We will begin our discussion with a model independent operator analysis, corresponding to very 
heavy mediation particles (such as a heavy Z' or squarks). In Section 2 we will introduce some 
representative four fermion operators supressed by a cutoff scale. We will then place limits on the 
strengths of these operators from the Tevatron mono-jet search. In Section 3 we will translate the 
Tevatron bounds to limits on direct detection cross section for different dark matter scenarios. In 
Section 4 we move on to introduce lighter mediators that are kinematically accesible at the Tevatron 
and find that these can either slightly enhance or severely weaken the Tevatron bounds. In Section 5 
we will discuss possible enhancements to the Tevatron dark matter search using the mono-jet p^ 
spectrum, and conclude. 

2 Operators and mono-jets 

Throughout this paper, we will assume a dark matter particle, x, as a Dirac fermion. The operators 
we will study are, 

Oz = -^^{qi,l,q){xri^x), 

Oa = ^J^ (9759) (X75X) , (3) 

Here we take q = u,d,s and turn on each operator one at a time (but results for a flavor universal 
operator will be easy to deduce), q'^ is the exchanged momentum and the suppression scale M is 
related to the mass of the particle whose exchange generates the four fermion operator. 

This is a representative set of operators that will generate a variety of dark matter scattering 
scenarios. Majorana dark matter will yield similar result (though for a Majorana spinor there are no 
vector interactions). Initially we will assume that the mediator is heavy and integrate it out, but in 
Section 4 we will discuss the effect of a light mediator. There are two additional operators x^^'^x P^v 
and WHxx appearing up to the dimension six level. While they are less constrained at the Tevatron, 
we leave their study and the study of operators involving the three heavy quark flavors to future 
work^ . 

Operator Oi leads to spin-independent coupling between the DM and a nucleus and can be thought 
of as arising from exchange of a scalar of mass M, O2 is similar but occurs through vector exchange. 



^ There is also considerable uncertainty in the heavy quark content of the nucleons for several of the operators under 
discussion [8]. 
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Figure 1: The constraints on the cutoffs of different operators from the CDF mono-jet search data at 
90% C.L.. 



Operator O3 is generated through axial-vector exchange and gives a spin-dependent couphng, and O^ 
could arise from exchange of a pseudo-scalar and gives a momentum dependent and spin-dependent 
DM coupling. Various combinations of these operators may be also generated by mediators charged 
under the SM such as squarks in super symmetry. 

2.1 Tevatron limits 

The CDF collaboration has performed a search for one jet events with large missing transverse en- 
ergy using 1 fb^^ of data [2]. CDF considered events with a leading jet pT and missing transverse 
energy both greater than 80 GeV. Events with a second jet with a pj- < 30 GeV were included but 
events with additional jets with transverse energy above 20 GeV were not. The number of observed 
events was 8449, a slight deficit compared to an expected background of 8663ib332. The standard 



model backgrounds are dominated by Z+jet, T^+jet with a missed lepton. QCD and "non-collision" 
background events contribute subdominantly to the background, but due to their high uncertainty 
they add a significant portion to the uncertainty of the background. The pT spectrum observed by 
CDF compares well with the expected background, however since the background uncertainty was 
only presented for the total number of events we will only use a simple counting experiment to place 
bounds. We will elaborate more on the potential of a multi-bin analysis in Section 5. 

We have generated signal events for the operators 01-4 using calcHEP [9] introducing a single 
operator at a time. We work with leading order parton level events, but have checked that including 
showering, hadronization and detector effects does not change our results significantly, see Section 5. 
Only the 80 GeV pT cut on the leading jet need be applied since we only generate leading order events 
and no additional jets are present. We add the signal events to the standard model expectation and 
require that the total number of signal plus background events do not disagree with the observation 
at more than the 90% confidence level. This allows us to place limits on the cutoff A = M/ ^g^gq 
which are shown in Figure 1. 

Note that the bounds on the various types of operators are remarkably similar, yet at low energies 
they will produce radically different direct detection signals. As expected, the Tevatron bounds are 
rather flat in energy for low dark matter mass, and start becoming weaker above 50 GeV, where 
kinematics and PDF suppression starts entering. Due to PDF's of different quarks, the up-quark 
operators always provide the most stringent bound. The constraints on the strange-quark operators 
are fairly weak, which justifies neglecting of even heavier quark operators. We note that, when the 
cutoff is close to the dark matter mass, the real constraints are subject to order one changes and we 
shall investigate this region further in Section 4. 

3 Direct detection limits 

In this section we translate the bound on the strength of the various operators into bounds on direct 
detection scattering rates. This enables us to plot Tevatron limits in the standard a — m^ plane. As 
we shall see, the Tevatron places competitive direct detection bounds for several interesting scenarios. 
We will summarize some of the interesting results here before presenting them in greater detail. In 
general collider bounds will be relatively more powerful in cases where direct detection scattering is 
suppressed, either by kinematics such as when dark matter is very light such as [10, 11, 12, 13, 14, 15] 
and the case of inelastic scattering [16], or by other suppressions such as when dark matter scattering 
is spin dependent or momentum dependent [17, 18]. This is because none of these suppressions will 
be effective for collisions at high energies. 



For instance the ~ 100 keV splitting between incoming and outgoing DM states in inelastic dark 
matter is not an impediment to their production at colliders. In addition, the Tevatron is not limited 
by the features of the DM velocity distribution in our galaxy and so is able to probe down to very 
low dark matter masses, well below the thresholds of direct detection experiments. Colliders will 
also be powerful in constraining direct detection scenarios in which a small fraction of dark matter is 
participating in scattering, but with an enhanced rate (such as [15] though in this particular model 
the interaction is mediated by a light boson which relaxes the bound). Additionally, vector couplings 
that are purely to strange or other sea quarks are not accessible to direct detection experiments and 
collider bound may be the only way to discover them experimentally. 

In our analysis we place explicit limits on single flavor operators, we comment on the bounds on 
combinations of operators at the end of this section. Flavor universal operators will have results that 
are close to the best single flavor operator (typically that of up). We note in passing that mono-jet 
searches will not be able to constrain models where the DM-SM coupling does not involve two neutral 
states from the dark sector. For instance resonant dark matter (rDM) [19] involves the DM state and 
a nearby charged state and so its Tevatron signal would instead be a jet, missing energy and a charged 
track. 

3.1 Spin independent 

The operators Oi and O2 induce a spin independent scattering of dark matter off of nuclei. To compute 
this scattering cross section off a nucleon, N = p,n, we will need to know the quark content of the 
nucleon {N\qT q\N) for each of these operators. At the nucleon level these operators become 

1 = ^1 A2 ' ^ ' 

,^Nq _ ,N {^rN){xi,x) 

^2 - Jq ^2 ' 

we consider the case M^ ^ q^ and A = M/ ^g^gq. The coefficients necessary to translate the quark 
level operators to the nucleon operators are given by [20, 21, 22] 



u 


= B2 


= 8.22 ±2.26, 


'd 


= B^ 


= 6.62 ±1.92, 


>P 


= B^ 


= 3.36 ±1.45 



(5) 

In extracting these conversion factors we have used the quark masses ratios m„/mrf = 0.553 ± 0.043, 
"^slf^d = 18.9 ± 0.8 and quark mass m^ = 105 ± 25 MeV[23]. 
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Figure 2: Left panel: the constraints on the spin-indepedent DM-proton scattering cross section. 
Relevant experimental bounds are shown as labeled. Right panel: the same as the left panel but for 
the constraints on the spin-indepedent DM-neutron scattering cross section. 



For the vector operator, 02, /« = /^ = 2 and /^ = /" = 1 and for all other quarks / = 0. Note 
this means that if the DM couples through vector couplings to second and third generation quarks 
only then it can never be discovered in direct detection experiments, but can be found using colliders. 
At low DM speed the leading contributions to the scattering cross section in each case are 



a 



Nq 



a. 



Nq 



7rA4 ^^5 ' 



(6) 

(7) 



where ^ is the reduced mass of the dark matter-nucleon system. The Tevatron limits on spin inde- 
pendent dark matter scattering for the various operators is shown in Figure 2. The recent results 
from CoGeNT [7], CDMS [3] and DAMA [6] with and without channeling are also shown in Figure 2. 
Note that the limits are slightly different for protons and neutrons simply because they are derived 
from proton rather than neutron collisions. The up-type and vector coupling operator are the most 
constrained operators. For dark matter with a mass below around 5 GeV, the mono-jet searches at 
CDF provide the world-best spin-independent bound. 

3.2 Spin dependent 

Models in which dark matter scattering is spin dependent are even more constrained by collider 
experiments. This is because SD scattering is suppressed relative to SI at low momentum transfer, 
because the scattering is not coherent over the whole nucleus, while there is no relative suppresion 
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Figure 3: Left panel: the constraints on the spin-dependent DM-proton scattering cross section for the 
up, down and strange (bottom to top sohd hnes) axial-vector operators. Relevant experimental bounds 
are also shown. Right panel: the same as the left panel but for the constraints on the spin-indepedent 
DM-neutron scattering cross section. 



between the two at high energies. Of the operators under consideration, spin dependent scattering is 
caused by the axial vector operator O3. For a complete list of all operators, see [24]. 

Again, in order to compute the DM scattering cross section off a nucleon, N = p,n, we will need 
{NpslN), leading to 



Nq 



with [21] 
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AJJ = 0.842 ±0.012, 



The total cross section is then 
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= K 


= -0.427 ±0.013, 


A? 


= A^ 


= -0.085 ±0.018. 
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The Tevatron limits on spin dependent dark matter scattering for the various operators are shown in 
Figure 3 along with limits from XENONIO [4], COUPP [25], PICASSO [26] and ZEPLIN III [27]. For 
the DM-proton spin-dependent scattering cross section (left panel) we have found that the Tevatron 
limits are stronger than any other direct detection experiments for all three operators. For the DM- 
neutron scattering in the right panel, the Tevatron limit is still the best for the up-type quark operator. 
Limits for a flavor universal operator are close to those of the pure up operator. 



3.3 Inelastic, Exothermic and streams 

We will now discuss constraints on several dark matter scenarios which have been proposed in the 

context of the DAMA modulation signal. The event rate in a direct detection experiment at a given 

recoil energy Er is proportional to 

dR r-- /(v) , 

-^^ oc n^aN / —;^dv , (10) 

where f{v) is the velocity distribution of dark matter and n-^ is the number density of the dark matter 
species in question. 

One interesting possibility is that if dark matter up-scatters from a ground state to a slightly 
excited one, the minimum velocity required for scattering is affected as 



' ^"^^^^+5 1 (11) 



2mTEji V /"T 

where rriT is the target nucleus mass, fiT is the target-dark matter reduced mass, and 6 is the mass 
splitting. This increase in the minimum velocity required to scatter causes the integral in equation (10) 
to be suppressed in iDM models (it also becomes more sensitive to velocity modulation, thus explaining 
DAMA well). In order to keep a certain event rate fixed, one would need to enhance the nucleon cross 
section by the corresponding factor. However, high energy collisions are not sensitive to the 0(100 keV) 
splittings between the dark matter states and thus collider bounds will have a relative advantage here. 

To explore the sensitivity of the mono-jet search to iDM we show the spin independent bounds 
again in Figure 4 now compared to the parameter space of iDM that explains DAMA at 99% C.L., 
taken from several studies. The DAMA preferred region depends sensitively on various assumptions 
to which our bounds are not sensitive. These range from astrophysical quantities such as the dark 
matter velocity distribution as well as experimental issues such as the DAMA energy resolution and 
the presence or absence of channeling in DAMA. In particular in Figure 4 we show the DAMA best 
fit region in the black contour region from Ref. [28, 29] at 99% C.L. for a fixed value 6 = 35 keV 
(for the lower mass region, corresponding to scattering off sodium) and 5 = 120 keV (for the higher 
mass region, corresponding to scattering off iodine). We also included another fit in the green dashed 
region from Ref. [30], where the parameter 5 has been treated as a floating parameter. All of these 
regions are also constrained by other direct detection experiments to various degrees, depending on 
the assumptions. We emphasize that the Tevatron is sensitive to a large part of the DAMA preferred 
region independent of these assumptions. 

Another class of models proposed to explain DAMA is exothermic dark matter [14, 15] which 
involves a light dark matter state which is de-excited upon collision with a nucleus. Here a larger 
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Figure 4: The constraints on iDM and exoDM with the DAMA best fit region for iDM as given in 
Ref. [30] (green, dotted), and in Ref. [28, 29] (black, solid). The DAMA best fit regions for exothermic 
DM taken from Ref. [14] is shown in red (small lower region) and from Ref. [15] in blue (small upper 
region), note that this model uses a light mediator, which will weaken the mono-jet constraints. 



cross section is required in order to enhance the modulation signal which is reduced in exothermic 
reactions. The allowed region from [14] lying below 10"^^ cm^ may also be probed in future Tevatron 
analyses. 

A third class of models in which collider bounds have a relative advantage are those in which 
direct detection signals are arising from a sub-dominant component of the dark matter halo. For 
example, one could exploit the fact that the number density of a thermal relic generically scales as the 
inverse of the annihilation cross section n oc l/{av) to argue that the rate count at a direct detection 
experiment, nav is a constant as the coupling strength of the thermal relic with matter is increased. 
As the relic couples more strongly it becomes less abundant while keeping the rate fixed. The Tevatron 
bound, which obviously does not depend on the number density of dark matter places a direct upper 
bound on a. For example, in [15] a light relic which accounts for about 1% of the dark matter energy 
density scatters with a cross section of a few x 10^'^® cm^, which is strongly constrained by the Tevatron 
mono-jet bound. 

In addition, some direct detection models for DAMA have relied on streams of dark matter that 
have a low velocity dispersion, but are sub-dominant contribution to the halo density. A stream model 
may also have an enhanced cross section, if the number density in the stream is below that of the 
generic halo. Examples of such possibilities that may be constrained by our bounds were explored 
in [29, 31]. 

It should be emphasized however, that all of these models may escape the Tevatron bound if the 
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interaction of the thermal rehc with matter is mediated by a hght state, as we shah see below. One 
may turn this statement around to conclude that if a dark matter model is directly detected in a 
region that violates the collider bound, then the dark sector is not simply a DM state, but contains a 
new light mediator through which the DM interacts with the standard model. 

In addition to considering the bounds on one operator at a time there are many models that 
predict more than one operator has non-zero coefficient. In these cases the direct detection cross 
section bounds presented above must be combined. If the direct detection cross section for a single 
quark flavor, q, is aq = dgji'^/Ag then in the case of multiple operators coupling quarks to DM the 
direct detection cross section is, 



Ei 



CTi 



4 Constraints on light mediators 

In placing the bounds in the previous sections we have imagined that the only accessible state from the 
dark sector is the DM itself, all other states associated with the dark sector are heavy [32]. However, 
for certain operators the cutoff scale, shown in Figure 1, is low enough to be probed at the Tevatron. 
In these situations it may be possible to produce the mediator that generates the four fermion operator 
directly. If the mediator couplings with the SM and the dark sector are weak, i.e. < 0(1), then the 
mediator mass is lower than the cutoff scale shown in Figure 1, further motivating consideration of 
mediators within the Tevatron's reach. Furthermore, recent cosmic ray excesses may be explained by 
a dark sector that contains a light mediator, M r^ 1 GeV, see for instance [33]. 

As discussed in the introduction, the ratio of the direct detection cross section to the mono-jet 
production cross section is proportional to 1/M^ when the mediator is light. So, as the mediator 
mass decreases, the constraints on the dark matter direct-detection from the mono-jet searches be- 
come weaker. For a sufficiently light weakly coupled mediator that satisfies the mono-jet bound 
the constraint on direct detection will not be competitive with those from direct detection experi- 
ments themselves. However, there is an interesting regime with both a light mediator and light DM, 
2m^ < M < s^'"^, where the production of DM in mono-jet events can take place through an on-shell 
mediator which in turn can decay to dark matter. In this situation the production of dark matter is 
a two body process rather than three body and so is enhanced by a phase space factor. Note that 
in this case the mediator could also have a substantial branching fraction to jets, leading to a di-jet 
invariant mass peak, though this is model dependent and will not be considered further here. 

For simplicity, we consider the mediator to be a SM singlet vector or scalar and consider the 

11 



5 10-38 

b 



10-*" 



U- 



'J\ GeV 



- — M;iOGeV 

Contact pp^^;_^-^ 
^-^GeV CoGeNT ^5^ 




0.5 1.0 5.0 10.0 

My (GeV) 



DAMA 
^(w. channeling) 

CDMS ~ 



50.0 100.0 



M . 1 n-33 - 



Q 10" 




0.5 1.0 5.0 10.0 

m^ (GeV) 



50.0 100.0 



Figure 5: Left panel: constraints on the spin- independent DM-neutron scattering cross sections for 
different mediator masses. Right paneh the same as the left panel but for the constraints on the 
spin-dependent DM-neutron scattering cross section. 



effects as its mass is lowered. In particular, we consider the case of both a 10 GeV and 100 GeV vector 
mediator, in both cases we consider the width to be 1% i.e. T = O.OIM. We leave the possibility 
of mediators that are charged under the SM, such as squarks, for future study (though their masses 
cannot be lowered below current direct bounds). 

As we alter the mass of the mediator we also alter its couplings to the SM and the dark sector, 
gq and gy^ respectively, so that the mono-jet production cross section satisfies the CDF bound. The 
effects of a light mediator for the case of SI DM, O2, and SD DM, O^, are shown in Figure 5. The 
weakening of the limits for light mediators is clearly seen, and the slight enhancement for the case 
where the mediator is produced on shell (M = 100 GeV and m-^ < 50 GeV) is also observable. 

4.1 Momentum dependent 

A particular model of dark matter that requires the introduction of light mediators is the case of 
DM that has a momentum dependent coupling to nucleons [17, 18]. These types of models provide a 
possible explanation for the DAMA modulation signal, but in order to do so require mediators less than 
10 GeV in mass. Although light from a collider perspective the masses considered are still sufficiently 
heavy that at direct detection experiments the mediator can be integrated out and an effective four 
fermion operator can be written. The axial-scalar operator O4 leads to momentum dependent and 
spin dependent dark matter scattering and at the nucleon level the operator is, 
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Figure 6: The constraints on the momentum and spin dependent model from mono-jet searches. The 
sohd hues are for a mediator with M = 10 GeV, while the dashed lines are a mediator with M = 1 GeV. 
The DAMA allowed region is shown in the green contours and is taken from Ref. [18]. 



where we have integrated out the mediator and A = Mj ^g^ q^ . In going from quark to nucleonic 
operators we introduce C^ = {N\qi'^^q\N) , which are determined in Ref. [20] (under the assumption 
that the large N^. limit is a reasonable approximation of QCD), 



C7P =141 ±34, C^ = -138±34, 
C^ =-137 ±33, C2 = 139±34, 
CP =-4.15 ±1.0, C^ = -0.94 ± 0.33 . 



(14) 



The differential cross section for DM scattering off a nucleon is given by 



da 



Nq 



K)\ 



(15) 



dcos9 327rA'^ (m^ ± mjv) 

where q is the exchanging momentum of the DM scattering off the nucleon. 

Following Ref. [18], we use a reference momentum, g^ef = 100 MeV, and compare the Tevatron 
constraints to the region of parameter space that best fits the DAMA result, taken from Figure 3(b) 
in [18]). The results are shown in Figure 6; we consider the cases of M = 1, 10 GeV. 

We see that the dilution of the Tevatron constraints by the light mediator means that momentum 
dependent dark matter with M = 1 GeV is not severely constrained by the mono-jet search. However, 
if instead the mediator is 10 GeV and has 0(1) couplings, then the lack of a mono-jet excess places 
strong constraints on the model and rules out the DAMA preferred region^, note that unlike previous 



This option may well be ruled out by other hmits. 



13 



ra^^lOGeV 



1000 



O 
o 



> 



10 



M=10GeV^ 

M=100GeV > 

Contact Op. ^ 




3.0 r 

2.5 ■ 
2.0 - 



o 

" 1.5 - 



C 

^ 1.0 

w 



0.5 
0.0 - 





A = 1 TeV 






"r''u XTfiX 1 A"* : 




- 




Calchep - 

Mad-partor 


' 


H 


Mad-pgs ^ 


■■t2 


- 


— T 


- 


tr.-- 


: 


^= 


-t^^Lz... ■ . ... 








100 



150 200 

Jet Pt (GeV) 



250 



300 



100 



150 200 

Jet Pt (GeV) 



250 



300 



(a) 



(b) 



Figure 7: (a) Comparisons of the shapes of the signal, the SM background and CDF measured events. 
The SM predictions are shown in the green and the CDF observed data are shown in red. (b) 
Comparisons of simulated signal events from two different Monte-Carlo tools and for the parton and 
the particle levels. The cutoff A = M/ ^g-^gq is chosen to be 1 TeV. 



cases, the constraints coming from the strange quarks are the most stringent. This is due to a small 
matrix element for the strange quark in equation (14). 

5 Discussions and conclusions 

It is worthwhile to consider possible improvements to the dark matter search at the Tevatron, and in 
the future at the LHC. Here we placed bounds on dark matter using only the total rate of mono-jet 
signal events above a certain px cut. An analysis that takes the spectrum shape into account may 
yield more powerful bounds. We show the spectral shape of the signal compared to the background 
in Figure 7(a). We find that the signal spectrum is somewhat harder than the background, especially 
when the messenger mass is much higher than the dark matter mass. We find that including showering, 
hadronization (using Pythia [34]) and a detector simulation (PCS [35]) does not change the signal 
shape significantly, particularly above 100 GeV, as is shown in Figure 7(b). This may allow us to 
place tighter constraints using a multi-bin analysis as compared with a simple counting experiment, 
since signal predicts more deviations in high px bins. However, this would require knowledge of the 
theoretical uncertainty on a bin-by-bin basis which is not presently available. Furthermore, a bound 
may be extracted from mono-photon events. 

In this work we show that the Tevatron mono-jet search places competitive bounds on dark matter- 



14 



nucleus cross sections relevant for direct detection experiments. In particular, the Tevatron limits are 
the current world-best for light dark matter, below a mass of 5 GeV. The Tevatron also sets the best 
limit spin dependent dark matter scattering. Various models built to explain the DAMA modulation 
signal such as inelastic and exothermic dark matter are also constrained by current Tevatron searches. 
In addition to considering dark matter that couples to quarks via contact interactions we have taken 
the possibility of light mediators, as motivated by cosmic ray excesses [33] into account. We find that 
the introduction of a light mediator of mass < 10 GeV alleviates the mono-jet bounds completely for 
most cases. This leads to an interesting conclusion - if a direct dark matter signal is established in 
a region that is in conflict with collider bounds, a new light state should be introduced to reconcile 
the data. The Tevatron is unable to search for DM above a few 100 GeV due to kinematics, an upper 
bound that will be raised at the LHC. The current powerful bounds using only 1 fb^^ motivates a 
dedicated analysis using more Tevatron data as well as future analyses at the LHC. 

Note added: While in the final stages of this work Ref. [36] appeared on the arXiv. They address 
similar issues, focussing on Majorana DM, but do not consider light mediators. 
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A Matrix Elements 

We present here analytic expressions for the mono-jet processes mediated by the operators in equation 
(3). We consider only the case of the mediator being sufficiently massive it may be integrated out and 
we take the initial state quarks as massless and label momenta as q{pi) + q{p2) -^ ffl?) +x(^i) +x(^2)- 
In the center of mass frame, the differential cross section for mono-jet production is 

daj 1 



dEidE2d cos edcpd-y 5Utt^E^, 



-IMI^ , (16) 



cm 
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where we have introduced final state energies and angles through, 

El 




, sm y cos 0, sm t/ sm 0, cos ( 



El — rn^ 



^2 = J E2 — m^ I — =,sin6'cos((^ + a),sm.9sm.[(j) + a), cos 6* I . (17) 



2-^X 



The matrix element for each operator is, 



2 _ 16 {kiM -'m^)p2-q 



,,.,2 32 (A;2.p2)(A;i.g) + (/ci.p2)(fc2-g) + "^^(P2.g) ,^„, 

l-^2p = T-TT , (19) 



9A4 pi.q 

32 (fc2.P2)(fei.g) + (fci.p2)i 
9A4 pi.q 



,,.,2 32 (fc2.p2)(A;i.g) + (A;i.p2)(A;2-g)-m2(p2.g) 

l-^sl = 7774 — , (^Oj 



2 16 (/ci.fc2+m^)p2.g 
'■^^' = 9A^ ^q • (^^) 

As before we define A = M/^g^Qg. 

B Relic abundance 

The DM abundance as measured by WMAP is Q^h^ = 0.1120 ± 0.0056. If the DM is a thermal relic 
its abundance after freeze-out can be related to its annihilation cross section [37] through, 

M' » '■"' " ;°° "1':"'', , (22) 

where the partial wave expansion of the annihilation cross section is ct' = a + hv"^ + . . . . For a typical 
weak scale WIMP 5* ~ 100 and xp = m^/Tp ~ 20 — 30 and the annihilation cross section should 
be (Tf '^ 3 X lO^^^cm^s^^. For each of the operators under discussion the annihilation cross section, 
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Figure 8: Abundance of DM if it were a thermal relic, the shaded region has abundance greater 
than the WMAP bound, fi-^/i^ = 0.1120. Reading from top left to bottom right the plots are for 
Oi, O2, O3, O4, and in each case we have considered only the coupling to one flavor of light quarks. 



expanded in relative velocity, v, is 
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We show the resulting relic abundance if the DM were a thermal relic in Fig. 8. It should be noted that 
the DM need not be a thermal relic, this is a model dependent assumption, and our model independent 
approach is sensitive to all regions. In particular, regions where the thermal relic abundance is below 
the WMAP bound are still viable regions to be probed by direct detection and collider searches. 
Furthermore, in the region where the abundance is above that of WMAP there may be additional 
annihilation modes, not useful for searches at the Tevatron, that would lower the abundance. 
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